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We study the influence of spin-orbit interactions on quasiparticle dispersions in two-dimensional 
electron and heavy-hole liquids in III-V semiconductors. To obtain closed-form analytical results, 
we restrict ourselves to spin-orbit interactions with isotropic spectrum and work within the screened 
Hartree-Fock approximation, valid in the high-density limit. For electrons having a linear-in- 
momentum Rashba (or, equivalently, Dresselhaus) spin-orbit interaction, we show that the screened 
Hartree-Fock approximation recovers known results based on the random-phase approximation and 
we extend those results to higher order in the spin-orbit coupling. While the well-studied case of 
electrons leads only to a weak modification of quasiparticle properties in the presence of the linear- 
in-momentum spin-orbit interaction, we find two important distinctions for hole systems (with a 
leading nonlinear- in-momentum spin-orbit interaction). First, the group velocities associated with 
the two hole-spin branches acquire a significant difference in the presence of spin-orbit interactions, 
allowing for the creation of spin-polarized wavepackets in zero magnetic field. Second, we find that 
the interplay of Coulomb and spin-orbit interactions is significantly more important for holes than 
for electrons and can be probed through the quasiparticle group velocities. These effects should 
be directly observable in magnetotransport, Raman scattering, and femtosecond-resolved Faraday 
rotation measurements. Our results are in agreement with a general argument on the velocities, 
which we formulate for an arbitrary choice of the spin-orbit coupling. 

PACS numbers: 71.10.-w, 71.70.Ej, 71.45.Gm, 73.61.Ey 



I. INTRODUCTION 

Semiconductor heterostructures offer the possibility of 
forming two-dimensional liquids with tunable sheet den- 
sity n s . In an idealized model where the carriers have 
parabolic dispersion with band mass m and interact with 
Coulomb forcesji the only relevant quantity is the di- 
mcnsionless Wigncr-Scitz radius r s = 1/ \JiTa 2 B n s , where 
as = h 2 e r /me 2 is the effective Bohr radius (e r is the di- 
electric constant). Since r s serves as the interaction pa- 
rameter, changing n s allows for a systematic study of the 
effects of the Coulomb interaction. In particular, proper- 
ties of quasiparticle excitations such as their dispersion 
and lifetime are significantly modified due to electron- 
electron interactions 

Great attention has been paid in recent years to band- 
structure effects involving the spin degree of freedom. 2 
The strength and form of spin-orbit interaction (SOI) 
can be controlled in two-dimensional liquids through the 
choice of materials, the type of carriers (electrons/holes), 
and details of the confinement potential. For example, it 
is possible to change the coupling constant with external 
gates i2r— In addition to detailed studies of single-particle 
properties, the problem of understanding the effects of 
SOI in the presence of Coulomb interactions is a topic of 
ongoing investigations. 

Quasiparticle properties in the presence of SOI 6 - - — 
have been examined primarily accounting for Rashb a 11 ' 12 
and/or Dresselhaus SOI ) 13 ' 14 which are dominant in elec- 
tronic systems. The SOI results in two distinct spin sub- 
bands, with two associated Fermi surfaces. The effects 
on quasiparticles are usually very small; at each of the 
two Fermi surfaces the quasiparticle dispersion 6 -! 7 -^ and 
lifetime 7 - - — are almost unaffected by SOI, except in the 



case of very large SOI coupling In fact, it was found 
with Rashba SOI that the corrections to these quanti- 
ties linear in the SOI coupling are absent i 7 - Although ex- 
plicit calculations are performed within perturbative ap- 
proximation schemes, notably the random-phase approx- 
imation (RPA),— ~— the SOI leading-order cancellation is 
valid non-perturbatively (to all orders in r s )r^ Similar 
arguments hold for other physical quantities! 16 ' 17 For ex- 
ample, values of the ground-state energy obtained with 
Monte Carlo simulations^ for up to r s = 20 could be 
reproduced with excellent accuracy by simply neglecting 
SOI corrections to the exchange-correlation energy— No- 
ticeable exceptions exhibiting larger SOI effects are spin- 
textured broken symmetry phases , 20 ' 21 non-analytic cor- 
rections to the spin susceptibilit y 22 ' 23 and the plasmon 
dispersion.— All these examples involve the presence of 
spin polarization (either directly22r— or indirectly^ -) , in 
which case the arguments of Ref. [l6| do not apply. 

Another interesting situation occurs when the SOI 
has a nonlinear dependence on momentum, thus cannot 
be written as a spin-dependent gauge potential . 17 ' 24 ' 25 
Then, the approximate cancellations mentioned above 
are not expected. Winkler has shown that the domi- 
nant SOI induced by heterostructure asymmetry is cubic 
in momentum for heavy holes in III-V semiconductors j 2 ^ 
This theoretical analysis was later shown to be in good 
agreement with magnctorcsistance experiments^ 7 - Re- 
cently, the relevance of this cubic-in-momentum model 
was supported by the anomalous sign and magnetic- 
field dependence of spin polarization in quantum point 
contacts i 28 ' 29 SOI quadratic in momentum can also be in- 
duced for heavy holes by an in-plane magnetic field . 30 ' 31 

That these band-structure effects can substantially 
modify standard many-body results is confirmed by re- 
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cent Shubnikov-dc Haas oscillation measurements in low- 
density hole systems ! 33 ' 34 For example, a surprisingly 
small Coulomb enhancement of the g-factor has been 
reported 3 - 3 - and puzzling results have also been obtained 
for the effective masses m± of the two hole-spin sub- 
bands (a = ±)M- A mechanism for the small g- factor 
enhancement was suggested in Ref. IU if SOI strongly 
distorts the ground state spin structure, the exchange en- 
ergy becomes ineffective in promoting full polarization of 
the hole system. 

In this paper, we focus on the effective masses m±. 
We note that m± are directly related to the quasiparti- 
cle group velocities v± at the Fermi surfaces. With this 
in mind, we find it more transparent to discuss the effects 
of SOI and electron-electron interactions on m± in terms 
of wavepacket motion, as illustrated in Fig.[TJ If an unpo- 
larized wavepacket is injected at the Fermi energy (with 
average momentum along a given direction), the sub- 
sequent motion is very different depending on whether 
the SOI is linear (n = 1) or non-linear (n = 2,3). In 
the former case (n = 1), corresponding to electrons, we 
have v+ ~ i>_ and the motion is essentially equivalent to 
the case without SOI. In contrast, for holes with strong 
SOI, v+ 7^ V- and the two spin components become spa- 
tially separated. The separation between the two spin 
components can become quite sizable, and it should be 
possible to observe such an effect with, e.g., Faraday- 
rotation imaging techniques . 35 ' 36 Electron-electron inter- 
actions, in addition to modifying the average velocity 
v = (f + + f-)/2 (an effect which is well-known without 
SO fe 37 ' 38 ), are also reflected on the velocity difference 
(i> + — V— ) between the two spin branches. 

With these motivations in mind, we pursue a study of 
the quasiparticle group velocity in the presence of linear 
Rashba and non-linear SOI. An accurate analytical treat- 
ment of the electron-electron interactions can be carried 
out at high density, and we restrict ourselves to this in- 
teresting limit. The regime of strong electron-electron 
interactions is much more difficult to treat (see Ref. [38| 
for a Monte Carlo study without SOI) but it represents 
a relevant topic for future investigations (r s ~ 6 — 12 in 
Refs. [H and [13). 

This paper is organized as follows: In Sec. [TT] we in- 
troduce a model Hamiltonian including a generalized 
SOI 3 ^ and we briefly review its non-interacting proper- 
ties, demonstrating that injected holes will separate into 
spin-polarized wavepackets. The Coulomb interaction is 
treated in Sec. IIIII by extending the classic treatment 
of Ref. H^. We describe several results of this screened 
Hartrec-Fock approach in detail, focusing on the quasi- 
particle group velocities and the interplay of spin-orbit 
and Coulomb interaction effects. A discussion of these re- 
sults is given in Sec.[lV] In Appendix 1X1 we give a general 
argument showing that corrections to the velocity from 
any linear-in-momcntum SOI can always be neglected to 
lowest order. Finally, a number of technical details are 
provided in Appendices [B] and [Cl 
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FIG. 1. (Color online) Motion of a wavepacket for a fixed 
time r. In (a) the SOI is not seen, because the difference 
in the velocities v± of the two spin components is too small. 
Panel (a) applies without SOI or to electrons with Rashba 
SOI (n = I). In (b) we show the effect of the non- linear 
SOI present in hole systems (n = 2,3). Since the two spin 
branches have significantly different velocities v±, there is an 
appreciable separation A = \v+ — V- \t. We also illustrate 
the propagation of non-interacting wavepackets (dashed) for 
the same time r. The effect of interactions (at high density) 
is to enhance both the average velocity (v > Vf) and the 
separation of spin-components: A > 5 = \v+ — v°_\t. 



II. NONINTERACTING PROBLEM 



In the high-density limit, r s is small and the system is 
well-described by a non-interacting single-particle Hamil- 
tonian. We consider here the following modelj 3 -^ includ- 
ing a generalized SOI with a linear-, quadratic-, and 
cubic-in-momentum dependence for n = 1,2,3: 



p* p»<r H 

Hq = h 17 

2 m 



(1) 



where p is the momentum operator, m the band mass, 
P± = Px i ipy, o± — &x i icrj,, and 7 the generalized 
spin-orbit coupling, with er the vector of Pauli matrices. 
The phy sical justification of this model has been given in 
Ref. [3f|: the n = 1 Hamiltonian contains the Rashba SOI 
present in electronic systems ; 11 ' 12 while n = 3 includes 
the analogous term generated by an asymmetric confine- 
ment potential for holesi 26 ' 27 ' 29 Finally, the n = 2 case is 
also relevant for holes, in the presence of an in-plane mag- 
netic fieldi 30 ' 31 While SOI terms of different form gener- 
ally coexist (see Appendix [A}, we assume here that one 
n value is dominant. This greatly simplifies the problem 
by preserving the isotropy of the electron liquid in the x- 
y plane and is a good approximation for several relevant 
situations. For example, it was found for hole o 26 ' 27 ' 29 
that the n = 3 term can be much larger than corrections 
to the SOI due to bulk inversion asymmetry^* 3 - 2 - Diago- 
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nalizing Hq yields the energy spectrum 

h 2 k 2 



2m 



(2) 



with a = ± labeling the two chiral spin branches. The 
corresponding cigenfunctions are^i 



D ikr 



(3) 



where k is a wavevector in the x-y plane, #k is the angle 
k makes with the x-axis, and L is the linear size of the 
system. 

It is useful at this point to introduce a dimensionless 
quantity g characterizing the strength of the spin-orbit 
coupling 
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-fh n k F 

Ep 



(4) 



where Ep = h 2 k F /2m is the Fermi energy without SOI, 
written in terms of the Fermi wavevector k p = ^J2ixn s = 
\f2/{aBr s ). While 7 has different physical dimensions 
for each value, n = 1,2,3, the dimensionless coupling 
g always gives the ratio of the spin-orbit energy to the 
kinetic energy. The coupling g thus plays a role analogous 
to that of r s for the Coulomb interaction. Taking 7 to be 
independent of the density (this is not always the case^) , 



then from Eq. ^ we have g 



(since 



Ep ot k F ). This suggests that in the high-density limit 
(r s — > 0), the effects of SOI are suppressed for electrons 
(n = 1), but remain constant (n = 2) or are enhanced 
(n = 3) for holes. This simple estimate already indicates 
a qualitative difference for holes relative to electrons. We 
will sec that this difference is indeed significant in the 
following sections. 

In the presence of SOI the two spin bands (a = ±) have 
different densities n±, giving the total 2D sheet density 
n s = n + + n_. Keeping n s (thus kp) fixed gives a con- 
straint 011 the Fermi wavevectors k± for the a = ± bands, 
fci + k 2 _ = 2k 2 F . We can then characterize the solution 
to this equation with a single parameter x- 



(5) 



The parameter x = ( n - — n + )/n s gives the chirality 
and is determined by both the SOI and electron-electron 
interaction^ We will assume for definiteness that 7 > 
such that x ^ and k + < For a fixed generalized 
spin-orbit coupling 7, the non-interacting value of x can 
be determined by setting the Fermi energies of the two 
bands equal, i.e., E+(k+) = E°_(k^). This equation im- 
mediately gives a relationship between g and Xi 
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(1 + x )n/a + (! _ x jn/a • 



(6) 



We denote the solution of Eq. ^ by Xo(d), which gives 
the non-interacting Fermi wavevectors k± = kp^Jl =p xo- 



Explicit expressions for xo (9) are given in Ref. [3ll We 
only cite here the small-g behavior, which is easily found 
from Eq. (f6j): 



Xo{g) 



(7) 



We arc mainly interested here in the properties of the 
quasiparticles and, in particular, their group velocities 
v± . The fact that the dispersion relation @ is a function 
only of the magnitude k is a consequence of the model 
being isotropic with respect to k, which allows us to dis- 
cuss the magnitude of the group velocity on the Fermi 
surfaces, 



1 dE° ± (k) 



dk 



hk° 
m 



\n-l 



(8) 



The above expression can be evaluated explicitly in terms 
of the non- interacting chirality Xo(.<?) : 



— =VlTXo(ff)± -g[l T Xo(g)] {n - 1)/2 
vp 2 

~l±|(n-l), 



(9) 



where vp = hkp/m is the Fermi velocity in the absence of 
SOI and in the second line we have expanded the result to 
lowest order in g, by making use of Eq. (|7|). Equation © 
shows that there is no relative difference in group veloc- 
ity for the a = ± bands when n = 1. In contrast, when 
n = 2,3 a sizable correction linear in g is present. This 
difference reflects itself on the evolution of an initially 
unpolarizcd wavepacket injected at the Fermi surface, 
schematically illustrated in Fig. [T] While for electrons 
the wavepacket remains unpolarized, for holes the two 
spin components spatially separate with time. Electron- 
electron interactions modify the non-interacting veloci- 
ties v±, but the qualitative difference introduced by SOI 
between electrons (n = 1) and holes (n = 2, 3) remains 
essentially unchanged. A similar behavior holds for other 
properties of the electron liquid as well>^ vanishing cor- 
rections to lowest order in g were found for the quasi- 
particle lifetime;^ the occupation^ and the exchange- 
correlation energy— if only Rashba SOI {n = 1) is in- 
cluded. 

Finally, another relevant quasiparticlc observable is the 
effective mass. This has recently been measured in hole 
systems through experiments on quantum oscillations.— 
This physical quantity is simply given by m± = hk±/v± 
and is thus essentially equivalent to v±. 



III. SCREENED HARTREE-FOCK 
APPROXIMATION 

Realistically, charged particles interact through the 
Coulomb potential so it is interesting to understand how 
the presence of SOI modifies the behavior of the quasi- 
particles. The fully interacting Hamiltonian is given by: 



2 trf e r \Ti 



(10) 
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where (for electron i) is as in Eq. (JT|) and the pres- 
ence in (|10p of a uniform neutralizing background is un- 
derstood. Although many sophisticated techniques exist 
to approach this problemfi^ the simplest approxima- 
tion to the quasiparticle self-energy is obtained by only 
including the exchange contribution: 

E„(k) = E° a (k) + K(k), (11) 

where 

re , n (1 + aa'cosn9') , 

K(k) = -}_^± — 2 J -n Wa ,V(\k-k'\). (12) 

k'a' 

Here, rik'± = 0(fc± — k') is the occupation at T = 
for the a = ± band, respectively, with Q(x) the Heavi- 
side step function. The first factor in the summation of 
Eq. (fl~2"j) . involving the angle 9' between k' and k, arises 
from the scalar product of the non-interacting spinors 
[Eq. ©], and takes into account the specific nature of 
the spin-orbit interaction (n = 1,2,3). 

To lowest order, V(q) is the Fourier transform of the 
bare Coulomb potential, 2ire 2 / (e r q). As is well known,— 
this form of the Coulomb interaction leads to an unphys- 
ical divergence in the quasiparticle velocity. By consider- 
ing an infinite resummation in perturbation theory, the 
screening of the Coulomb interaction removes the diver- 
gence, e.g., in the RPA approximation. Finally, by ap- 
proximating the dielectric function in the effective inter- 
action by its zero-frequency long- wavelength limit, the 
RPA self-energy gives Eq. (TT2|) with 

V {q) = -Z%*L-. (13) 
q + V2r s k F 

This screened Hartree-Fock approximation with SOI,— ^ 
notwithstanding its simplicity, becomes accurate in the 
high-density limit. 

A. Renormalized occupation 



where we have rescaled the wavevectors k = kpy and 
defined y± = \/T^f~x~. The integral on the right-hand 
side is the correction from the exchange term and we 
have verified that Eq. ([15)) gives Eq. © for r s =0. We 
note that, for a given value of the SOI, x( r siff) enters in 
a rather complicated way in Eq. (|15p. being involved in 
the integration limits of the exchange term as well as the 
integrand. In practice, instead of solving Eq. (|15[) for x, 
it is convenient to evaluate g for a given value of x an< i 
numerically invert the function g(r s ,x)- 

The values of k± = kpy± were obtained in Ref. l3ll 
through a different procedure, i.e., by minimizing the to- 
tal energy (including the exchange contribution) of non- 
interacting states42 Although both methods are unreli- 
able at r s > 1, they both become accurate in the high- 
density limit, r s < 1. In fact, the only difference in the 
two approaches is due to the presence of the Thomas- 
Fermi screening wavevector in Eq. (|15j) and, by neglecting 
\J 7 ir s in the denominator of the second line, the equation 
from the variational treatment is recovered. In particu- 
lar, expanding Eq. (|15l) at small r s and g gives the same 
result found in Ref.l31t 

X^.^l-ifig^J. (16) 

Details of the derivation of Eq. (fTB")) are given in Ap- 
pendix [Bj Two salient features of Eq. (fTS)) are>21 (i) For 
n = 1 there is no correction to the noninteracting result 
x{ r s,g) — .<?• On the other hand, the linear dependence 
on g is actually modified by electron-electron interactions 
at n = 2, 3. (ii) The effect of the electron-electron inter- 
actions is a reduction of x( r s, g) from the non- interacting 
value. This result could be rather surprising, having in 
mind the well-known enhancement of spin polarization 
caused by the exchange energy^ (when the spin-splitting 
is generated by a magnetic field). However, x does not 
correspond here to a real spin-polarization, which is zero. 
Instead, x is simply related to the population difference 
of the two chiral spin subbands. 



By including the SOI, we can verify that Eq. (fT2"|) gives 
the correct high-density behavior for the Fermi wavevec- 
tors k±. These are modified by electron-electron inter- 
actions from their non- interacting values fc^i 31 ' 41 From 
Eq. (|12[) . k± can be obtained by equating the chemical 
potentials in the two spin branches: 



E+(k+) = E_{k.) 



(14) 



After taking the continuum limit, this equation is rewrit- 
ten in dimensionless form as follows: 



V2^Jo 2tt7 
y(a + a 1 cos n 



^/y 2 + yj - 2yy a cos 9 + \/2r, 



(15) 



B. Quasiparticle velocity 

In the screened Hartree-Fock approximation, the group 
velocities at the Fermi surfaces are given by: 



v± 



ldE k± 



h dk 



k—k± 



!^±±n 1 {hk ± ) n - 1 
m 

(1± cr' cos n9') 
- L 2^2 



dk 



nik-k'D 



(17) 

k— k± 



In general, we can discuss all corrections to v± by intro- 
ducing the following notation: 



v± 

VF 



1 + Sv(r s ) + Sv° ± (g) + Sv ± (r s ,g), (18) 
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where 5v(r s ) is the (spin-independent) correction due to 
electron-electron interactions at g = 0, which has been 
the subject of many theoretical and experimental studies 
(see, e.g., Ref. U H3, I37M391 and references therein). In 
the approximation (fl7|) . it is given by^- 



V2r s r 2 s r s {l-r 2 s )cos\x- l {V2/r s ) 



V2tt y/1 -r 2 s /2 



(19) 



As is known ji this approximation gives the correct lead- 
ing behavior at small r s : Sv ~ — (r s lnr s )/(-\/2vr). The 
second nontrivial term in Eq. (|18[) is the non-interacting 
correction purely due to SOI 



Sv° ± (g) 



Up 



1. 



(20) 



which only depends on <? [see Eq. (JSJ)]. Finally, <fo±(r s , p) 
collects all remaining corrections. 

A pictorial representation of the physical meaning of 
the three terms is shown in Fig. [1] for an unpolarizcd 
wavcpacket injected at the Fermi energy. At high den- 
sity, as seen in Eq. (fT9f and illustrated in Fig. QJa), the 
group velocity is larger than without electron-electron 
interactions. In Fig. [T{b) we depict the generic situ- 
ation with SOI. In the presence of SOI, the two spin 
branches have different group velocities. An initially un- 
polarized wavepacket then splits into its two spin com- 
ponents. Both the SOI and electron-electron interactions 
influence the relative velocity (v + — The separation 
after a time r is 6 = |in_ — v_\r for the non- interacting 
case and is modified by Sv± with electron-electron inter- 
actions: A = \v + — v_\t [see Fig.QJb)]. 

The presence of 'interference' terms, Sv±(r s ,g), be- 
comes clear from Eq. (TIT)) . These terms are due to the 
interplay of many-body interactions with SOI. A first 
contribution to Sv±(r s , g), which we refer to as the 'self- 
energy contribution', comes directly from the exchange 
integral [third line of Eq. (|17|) ]: due to the presence of 
two distinct Fermi wavevectors k±, the result obviously 
contains corrections to Eq. (fT9| which depend on g (in ad- 
dition to r s ). A second contribution to 5v±(r s ,g) comes 
indirectly from the non-interacting part and we refer to 
it as the 'rcpopulation contribution'. Since the Fermi 
wavevectors k± are modified from the non-interacting 
values fc°, the second line of Eq. (|T7f gives a result 
distinct from v±. The sign of this repopulation con- 
tribution is easily found by noting that, as discussed in 
Sec. lIII Al the exchange energy reduces the value of \ (f° r 
n = 2,3). This corresponds to an increase (decrease) of 
k + (k-), i.e., a positive (negative) correction to v + (i>_). 
The effect would be to enhance the difference in veloc- 
ity between the two branches, as illustrated in Fig. Qlb) 
(A > 5). However, to establish the ultimate form and 
sign of 5v± (r s , g) requires a detailed calculation of both 
self-energy and repopulation contributions, which is pre- 
sented below for some interesting cases. The total spin- 
dependent part of the velocity, 8v*± + 6v± , can then be 
compared to the simple non- interacting effect, Sv±. 



C. n — 1: higher-order corrections in g 

We begin the analysis of Eq. (fT7| by rewriting it in a 
more explicit way. To this end, we use 



^nik-k'i) = -k~v(ik. 



dk' 



k'l 



(21) 



where k = k/k. This allows us to integrate Eq. (fT7|) by 
parts, which leads to: 



v ± ,9 n-1 

— = V± ± »ny± 
vp 2 



271 dff_ 

8^ 



f 27r d9' 



\/2cos0'(l ± a cos n6')y cr 
y 2 ± +yl- 2y±y <T cos 6' + r s y/2 
y/2n sin nO 1 sin 9' dy' 
y\ + y' 2 - 2y'y± cos 9' + r s \/2 



— Lt\ + L2 + £3 



(22) 



Notice that, in the integration by parts of Eq. (|T7|) in the 
continuum limit, two types of terms enter: those corre- 
sponding to the second line of Eq. (|22|) (£2) involve the 
derivative of rtk'±- This results in a delta function in the 
dk' integral, which can then be easily evaluated. Thus, 
only the integral in d9' is left. The second type of term 
involves the derivative of cos n9': 



a d Ti 
k • — — cos n9' = — sin n9' sin 1 
ok' k' 



(23) 



and indeed this angular factor appears in the third line 
of Eq. (ED) (L a ). 

Eq. ([22"]) can always be evaluated numerically, after ob- 



taining the values of \ (thus y± = ^/l =j= X) from Eq. (fT5|) . 
By specializing to the small g, r s limit for n = 1 SOI, it is 
known that the linear-in-g correction to v± vanishesii^ 
Thus, an expansion to second order in g has to be per- 
formed, which has been done in Ref. @ in the context of 
the RPA treatment of the quasiparticlc properties. To 
verify the validity of the simpler screened Hartree-Fock 
procedure, it is interesting to perform the same expansion 
for our Eq. (|2^|) . 

In fact, the two calculations bear some similarities 



since Li is the same as the boundary term B } 



boundary 



of 



the RPA treatment, sec Eq. (69) of Ref. Thus, we can 
borrow the expansion in small g,r s : 

(24) 

where terms of order 0(r 2 ,r s g 2 ) have been omitted. To 
the same order of approximation, we have for L\\ 
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i-ir; (" = 1), 



(25) 



which can also be obtained by expanding Eq. © (with 
n = 1). Since x only receives 0(r s g 3 In g) corrections 
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from electron-electron interactions^ it is sufficient to use 
the non- interacting value, Xo(s)j to this order of approx- 
imation. Thus, the rcpopulation induced by electron- 
electron interactions has a negligible effect on v± in this 
case. The situation will be different for n = 2,3. Ex- 
panding L3 for small r s , g yields 



r ( n 

L 3 ^±-^—g; ("- = 1), 



(26) 



which cancels the linear-in-g term of Eq. (|2"4")l . as expected 
for n = 1. 

We note that in the screened Hartree-Fock approxi- 
mation we are able to obtain an analytic result for the 
leading term of L 3l and some details of the derivation can 
be found in Appendix O In contrast, in the RPA treat- 
ment of Ref. [D it was not possible to expand the more 
involved corresponding term, i?i n t, in a fully analytic 
fashion. The cancellation of the linear-in-g contribution 
was indicated by a general argumcntii^ and confirmed 
through numerical study. Additionally, the absence of 
higher-order terms which modify Eq. (|24p was inferred 
numerically in Ref. @. Although the final results of both 
approaches (RPA and screened Hartree-Fock) agree, the 
situation is clearly more satisfactory within the frame- 
work of Eq. ([22]) . since cancellation of the linear term in 
Eq. (|2"4"|) can be checked exactly and the expansion can 
be carried out to higher order systematically. The final 
result, computed to higher order in g, reads (for n = 1): 



v± r 

V F 



l + 6v(r s ) 



y/2r 



I6tt 



Sv° ± (g) 



\n g -+ 3 - 

8 2 



6 



319 
~20~ 



(27) 



The second line of Eq. (f2"T)) represents the expansion of 
Sv± (r s , g) including all terms up to 0{r s g 3 ). For n = I, 
the non-interacting result gives the same velocity for both 
spin branches [8v±(g) is actually independent of ± for 
this particular model]. The 0(r s g 3 ) term in Eq. (f2"T|) is 
therefore quite interesting. It shows that a small differ- 
ence in velocity exists. This is a genuine effect of electron- 
electron interactions. 

The accuracy of Eq. (|2"7) can be seen in Fig. [5] it 
becomes accurate at very small values of r s (as shown 
in the inset) while at larger more realistic values of r s 
it still gives the correct magnitude of the effect. The 
numerical example of Fig. [5] also shows that Sv±(r s ,g) 
is very small. In fact, it is generally much smaller than 
the non-interacting correction: since 5v±(g) ~ —g 2 /8 
[see Eq. (|25[) ] the r s g 2 hig term becomes larger only if 
r s In g 3> 1. This condition is only satisfied at extremely 
small values of g (if r s is small as well) at which SOI 
effects are hardly of any relevance. Thus, Eq. (|2"T|) implies 
that the effects of SOI and electron-electron interactions 
are essentially decoupled for n — 1 . This picture changes 
substantially for hole systems {n = 2,3), as we show in 
the next section. 
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FIG. 2. (Color online) Thick solid lines: the corrections 
Sv±(r s ,g) evaluated by numerical integration of Eq. (|22|l . 
Dashed lines: approximation to Sv±(r a , g), given by the sec- 
ond line of of Eq. (|27[) . In the main plot, we take r s = 0.1, for 
which Eq. (|27[) is not very accurate. Very good agreement is 
obtained at small values of r a (inset, with r s = 0.001). The 
thinner solid line in the main plot is the total correction due 
to spin-orbit coupling: Sv±(g) + Sv±(r s , g). It is dominated 
by the non-interacting effect and the spin-splitting is not vis- 
ible. The value of Sv±(g) at g = 0.1, outside the plot range, 
is ~ 1% (all corrections are in units of vf)- 



D. Corrections to the velocity for n — 2, 3 

Wc now apply the discussion of the previous section to 
the SOI more appropriate for holes and point out some 
important differences. Expansions for small r s and g are 
given in Appendix [C] From Eqs. (|C16p and (|C17[) and 
using x = 9 + 0(r s g, g 3 ) [Eq. (fTB]) ]. we find the following 
expressions for the self-energy contribution. For n = 2: 



L 2 + L 3 ~ Sv(r s ) 
and for n = 3: 



V2r s 



4tt 



± 9 - 
3 



(28) 



L 2 + L 3 



5v{r s 



4-7T 



± — g 
15 y 



U91n£ 



613 s 

4 V"~8 ' lo" 

(29) 

At variance with the case of n = 1, the linear-in-<7 term 
does not vanish here. Thus, we find an appreciable cor- 
rection to the velocity. This correction has opposite sign 
in the two branches and is positive for the + (higher- 
energy) branch. 

A second contribution to the g-linear correction comes 
from L\. Expanding L\ in terms of x and using Eq. (|16[) 
gives, for n — 2, 



L 1 ~l + 6v° ± (g) + 



and for 11 = 3: 



U ~l + 6v° ± {g) + 



V2r, 

47T 



V2r t 
Air 



2 

±-g 

3 y 



1G 



32 , 

± 15- 9+ 15 5 ' 



(n = 2), (30) 
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FIG. 3. (Color online) Plot of 8v±(r s ,g) for n = 3 (solid 
curves) and n = 1 (dashed curves) as a function of the SOI 
strength g. We have taken r s = 0.3 in both cases. ± indicate 
the spin branch of each curve. 



Thus, the repopulation contribution is present in this 
case and has the sign discussed at the end of Sec. MI Bl 
(it is positive for the + branch). 

As it turns out, the self-energy, repopulation, and non- 
interacting contributions to the velocity have the same 
sign. The three contributions therefore have a cooper- 
ative effect in enhancing the difference in velocity be- 
tween the two spin branches. Of course, based on the 
high-density theory presented here, we cannot tell if this 
conclusion holds at all densities. We also note that the g- 
linear term of the self-energy correction [Eq. (|28|) or ([29]) ] 
is always half of the corresponding repopulation correc- 
tion [Eq. (f30j) or ([3ljl], Again, we have not investigated 
if this curious relation only occurs within this approxi- 
mation scheme or if it is more general. 

Finally, we give the complete result for n = 2 



VF 



1 + Sv(r s ) + Sv° ± {g) 



V2r 

47T 

and for n = 3 
v± 



±g + g 2 (In! + 3 



(n = 2), (32) 



Vf 



I + Sv(r s ) + Svl(g) 



47T 



, 8 2 ( n , 9 869 

±-g + g 9 In - H 

5 y y \ 8 120 



(n = 3), (33) 



and show two numerical examples in Figs. [3] and [4] Fig. 
[3] is a comparison of Sv±(r s ,g) for n = 1 and n = 3: it 
is clear that the dependence on g is very weak for the 
electron case (n = 1) and the magnitude of 6v±(r s ,g) 
is much larger for holes (n = 3) with SOI of comparable 
strength. Fig. 2]shows the separation A between the two 
spin components of an initially unpolarized wave packet 
for a fixed travel distance of 1 /iin: 



A 



V-\ 



v + + V- 



x (1 /im). 



(34) 



FIG. 4. (Color online) Thick solid line: separation A be- 
tween the two spin components of a hole wavepacket (n — 3) 
at a fixed drift distance of 1 /im, see Eq. (|34|) . The dashed 
line is the non-interacting value Aq. The thin solid line is 
n = 1, indistinguishable from A = 0. We have assumed here 
r s — 0.3. The inset schematically illustrates the definitions of 
A and Aq. 



In addition to including both n = 1 and n = 3, in 
Fig. @] we plot the non-interacting value Ao, obtained 
by substituting v± — > v± in Eq. (|3~4")l . For n — 1 the 
non- interacting velocities v± are the same (Ao = 0). 
The effect of electron-electron interactions is not visi- 
ble. Thus, the wavepacket remains essentially unsplit 
and unpolarized (A ~ 0) and the only significant influ- 
ence is on the average velocity (w++w_)/2, from Eq. (fT9|) 
[see Fig. Ufa)]. In contrast, a large splitting is found 
for n = 3 where A, Ao can reach a large fraction of the 
traveling distance. Since spin-polarized ballistic trans- 
port is observed in two-dimensional hole systems on /xm- 
scales (e.g., in spin- focusing experiment a 28 ' 29 ) and the 
typical wavepacket traveling time is 1 — 50 ps (depend- 
ing on the density), Fig. 2] suggests that a direct op- 
tical imaging of the wavepacket separation should be 
within reach of femtosecond-resolved Faraday rotation 
measurements i 35 ' 36 

As for electron-electron interaction effects, we see in 
Fig. |4] that a visible difference between A and Ao exists. 
The difference is quite small, due to the fact that we are 
considering here the weak-coupling limit, and all the in- 
teraction corrections are proportional to r s < 1. Ao is 
modified here by ~ 2 — 3% and we collect some represen- 
tative values in Table [H It can be seen in Table [J that, 
while Ao does not change with r s , the interaction correc- 
tion A — A grows at lower densities (see Table |T|. Since 
experiments on hole systems can reach values as large as 
r s = 6 — 12j 33 ' 34 it is reasonable to expect significant ef- 
fects from Sv± in this low-density regime. As a reference, 
in electron systems, Sv changes from ~ 5% to —30% for 
r s ~ 1 to 

It might be surprising to see A < Ao in Fig. [4j This is 
due to the interplay of two competing interaction effects. 
It is true that, as discussed already, the difference be- 
tween v± is enhanced by Sv±. The two spin components 
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SOI g r s Ap (nm) A - Ao (nm) 

n=l 0.05 0.5 0.001 

n=3 0.1 0.1 203 -4.7 

n=3 0.1 0.3 203 -6.2 



TABLE I. Separation of wavepackets after 1 /im drift dis- 
tance for electrons (n = 1) and holes (n = 3). The non- 
interacting value and the electron-electron interaction correc- 
tion are listed, see Eq. (|34[) and Fig. [3] Typical values for 
g (~ x) are li, 0-05 for electrons^ and < 0.2 for holes.— ^ 
We have assumed in this table that typical values of g are 
independent of r s and have used high-density values r 3 < 1. 



therefore split faster in the interacting case. However, 
the situation shown in Fig. U is distinct from that shown 
in Fig. Q] since a constant traveling distance (and not 
time r) is assumed. At high density the mean velocity 
(i>+ + v-)/2 is greater for the interacting gas (Sv > 0), 
which allows the wavepackets less time to separate. As 
it turns out, the latter effect is dominant in Fig. 

Interestingly, the sign of Sv changes at low density^ 
This would imply a cooperative effect of Sv± and Sv on 
A if Sv± does not change sign. Unfortunately, to infer 
the behavior of Sv± at low density requires a much more 
sophisticated approach. 



cle dispersion. In addition to being relevant for trans- 
port measurements of the effective mass, the significant 
difference in group velocity at the Fermi surface of the 
two spin branches could also be addressed by Raman 
scattering ex per iments, demonstrated for electron sys- 
tems in Ref. |43|, or via time-resolved Faraday-rotation 
detection of the spin-polarizatio m 35 i 36 A similar discus- 
sion should hold with n = 2,3 for other physical obscrv- 
ables and many-body effects. For example, studying the 
compressibility^^ in the presence of SOI and extend- 
ing the n = 1 discussion of the quasiparticle lifetime" 
would also be topics of interest. Finally, the problem of 
including in our framework a more general form of SOI 
than Eq. (|T|) is clearly of practical relevance. However, as 
discussed in Appendix \^ we expect that our qualitative 
picture on the different role of linear and non-linear SOI 
remains valid also in this more general situation. 
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Appendix A: General spin-orbit coupling 



IV. CONCLUSION 

We have presented a theory of the quasiparticle group 
velocity at high density, in the presence of SOI of dif- 
ferent types. Contrasting the behavior of electron and 
hole systems, we find several intriguing differences. Wc 
have shown explicitly that the lowest-order cancellations 
of SOI effects occur only for the electronic case, when 
the SOI is approximately linear in momentum (e.g., a 
strong Rashba or Dresslhaus SOI is present). On the 
other hand, SOI terms non-linear in p are often domi- 
nant in hole systems. 26 ' 27 ' 29 Thus, larger effects of the 
SOI and a non-trivial interplay with electron-electron in- 
teractions are expected for holes on general grounds. 

As an important motivation for future theoretical stud- 
ies, hole liquids can be realized in the laboratory with 
strong SOI and large values of r s . For example, the spin- 
subband population difference at zero field is of order 
15 - 20% in Rcf. H3, with n s ~2-4x 10 10 cm- 2 . With 
a hole effective mass m ~ 0.2too in GaAs these densities 
correspond to r s ~ 9 — 12. For electrons, materials with 
strong SOI typically have small effective masses, which 
results in much lower values of r s . A diluted electron liq- 
uid with n s ~ 2 x 10 10 cm -2 gives r s ~ 1.2, using InAs 
parameters (to = 0.023mo). 

Discussing the large-r s regime of holes would require 
extending many-body perturbation theory^— or Monte 
Carlos approaches, so far only applied to linear SOI. 
The high-density regime studied here would represent a 
well-controlled limit of these theories for the quasiparti- 



A general SOI contains terms with linear and non- 
linear dependence on momentum, and does not neces- 
sarily have the isotropic form assumed in Eq. ([1]). For 
definiteness, we suppose that there is no magnetic field, 
so that quadratic terms are not present: 

H = Y.§^ + H T + HT + H cX _ eh (Al) 

i 

where Hf°, and -ff c i-ci give, respectively, the linear- 
in-momentum spin-orbit, the cubic-in-momentum spin- 
orbit, and electron-electron interactions. To show that 
Hf° always has a small effect, we consider the uni- 
tary transformation H' = e~ He with S defined by 
[S,£,p?/2m] =Hf, giving 

H' ~ |- + Hl° + ff cl _ cl - [S, Hf> + flf] . (A2) 

In the same transformed frame, the velocity operator of 
electron i is given by v£ = e~ (dH/dpi)e s , which yields 

, pi dm° r „ dm° ... 
v ~^ + — JL_ 5 i + J_ . A3 
to dpi dpi dpi 

In deriving Eq. (|A2j) we have used the fact that 
[S, flei-ei] = 0. For example, for an isotropic SOI as 
in Eq. ©, 

S = i— y^X x i a v,i - yi a x,i)- (A4) 
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The property [S, i? e i-ei] = is valid for a general 
linear-in-momentum SOI, including a combination of 
Rashba and Drcssclhaus SOL— However, the same iden- 
tity [S, H e \-ei] = does not hold for a transformation 
with [S, J2iPi /2ti] = Hf, i.e., a transformation that 
is designed to remove the non-linear component from 
the non- interacting Hamiltonian. In writing Eq. (|A3|) . 
we have used [5, Pi/m] = dHf/dpi, which implies a 
cancellation of the Hf contribution to v' { to lowest or- 
der. Again, this cancellation is only valid for linear-in- 
momentum SOI. 

By introducing dimensionlcss couplings associated 
with Hf° 3 , in direct analogy with Eq. we see that the 
commutators in Eqs. (|A2[) and (|A3|) are of quadratic or 
bilinear order in the couplings (~ g\ and ~ gigs)- This 
indicates on general grounds that has the largest ef- 
fect on the quasiparticlc velocity if g\ < g 3 <C 1. In 
this case, if we are interested in lowest-order effects, 
we can neglect both anticommutators in Eqs. (|A2[) and 
(|A3[) , which is equivalent to neglecting Hf° in the origi- 
nal Hamiltonian (|A1[) . Thus, to leading order all results 
we report for the quasiparticle velocities due to a pure 
cubic-in-momentum spin-orbit interaction also apply in 
the case of a mixed linear-plus-cubic spin-orbit interac- 
tion (with the caveat that we consider only the isotropic 
form of cubic SOI) . 



Appendix B: Derivation of Eq. (Il(>|) 

As discussed in the text, at high density we can neglect 
the y/2r s in the integrand of Eq. (fT5|) (second line). For 
small g, the value of x is also small and we can perform an 
expansion of the exchange integral. First notice that the 
constant term at x — is missing, because the integration 
limits simply become y± = 1 and the integrand vanishes 
upon the summation on a, a'. Therefore we only need to 
compute the linear term in x : 



C~* 7 — f V °' ( cr + cr ' cos »% ( fy 

dX ^ Jo 27T J ^y2+yl-2yy a COS 9 



d6 
2^ 



2(1 — y cos 0)ydy cosnO 



o (1 + y 2 — 2ycos9) 2 



(Bl) 



and after evaluating the dy integral in the square paren- 
thesis, Eq. (|Bip gives 



r 2lT de_ 

Jo 2^ 



2 cos 9 — cos 710 — 1 



-2 In 1 



sin ■ 



cos 9 



7T ^ 

J=0 



2.? - 1 



(B2) 



Equation ([16]) is then easily obtained from Eq. (| 15[) by 
neglecting all the cubic terms in the small parameters 
5:Xi r s [ e -S-i the left side of Eq. (TT5]) is {y\_ + y 2 _)g = 
2g + 0(g X 2 )]. 



Appendix C: Small r s ,g expansions 

We give in this appendix some details on the expan- 
sions of Eq. 



V -^=y ± ±lny^ + ^(I 1+ I 2+ I 3 ), (CI) 

VF 2 lD7T 

where we have split L2 in its two contributions (ii 2 refer 
to a = ±) and I3 corresponds to L3: 



f 2n V2y±cose(l + cos nb 
1\ = I d9- 
Jo 



V2y±sm9/2 + r s 

f 27T , n V2y T cos 6»(1 - cos n9) 
1 2 = / dO — , 

Jo v 1 _ y+y- cos » + r « 

f 27T , n f V ~ , 2nsinn9sm0 
h = ± / a0 / ay , =- 

"'° Jy+ Jy± + y 2 - 2yy± cos + r s V2 



(C2) 
(C3) 



Notice that these integrals have an explicit dependence 
on r s and y± = ^/l =p So, it is easier to perform first 
the expansion in the two small parameters r s ,x- The 
final results in the main text are given in terms of the 
physical couplings of the hamiltonian: r s and g. Those 
final expression are easily obtained by substituting the 
value of x i n terms of r s and g (x — g in first approxi- 
mation). 

The first integral, Eq. (|C2[) . can be evaluated exactly. 
In particular for n = 1 we obtain 



— + 8^(1 -5 2 ) 



fl - 3S 2 + 25 



4 ^ tanh 1 a/1 - 5 2 



(C5) 



where S = r s /\/2y±. This expression can then be easily 
expanded in r s , x and an analogous procedure is followed 
for n = 2,3. To lowest-order in r s , we can set r s = in 
I2 and I3. Similarly to the I\ angular integral above, 
the dO integrals of I2 and ^3 at r s = can be computed 
analytically for n — 1,2, 3. For I2 this yields directly the 
desired function of x- F° r -^3 we still need to perform a 
last integration in dy. Since the integration region is of 
size ~ x around y = 1, wc can expand the integrand in 
the small parameter (y— 1) and perform the integration in 
dy order- by-order, which allows us to extract the leading 
terms of the expansion in x- For n = 1 all this gives 



5h ~T4x-2 X 2 T3X 3 , 



SI 2 
SI 3 



4 

, 16 4 2 , 3 





13> 


4 


(-1- 


-1) 




" 6 ; 









In 



389 \ 



3 16 120/ 



(C6) 
(C7) 
(C8) 



where only the corrections SI a = I a (x) — I a (x — 0) are 
listed, since terms independent on x simply give the small 
r s expansion of the well known Eq. (I19[) . Here, terms of 
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order 0(r s x,X 4 ) are omitted, while it is interesting to The final results are for n = 2: 
keep the 0(x 3 ) terms, since they give the leading spin 
splitting. Indeed, it is easily checked that the linear terms 
cancel 

|>w(| + i n §)4(i4 + ^).<c9) £ 4 ^ +4x , (ta | +2)i (C16) 

a— 1 

which immediately leads to Eq. (|2"7|) . 

For n = 2, 3 we can proceed in a similar way. The spin 
splitting appears now already to linear order in x- By 
keeping the first subleading correction in x we have for 
n = 2: anc j f or n = 3 : 

«!~T4x-2x 2 , (CIO) 

„ s , ± | z + ^, (C12) ^^^(.kj^). (C17) 

and for n = 3: 



Q = l 



«!~ T 4x-2x 2 , (C13) 

(5/2 ^±^ x + x 2^ 91n X + ^\ (cw) 

105 A A V 8 42 / ' V ; 

, fi From Eqs. (|C16|) and (|C17[) we immediately obtain 

SI 3 ~ ±— X + ^X 2 - (C15) Eqs. m and ®. 
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